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There have been enough successful experiments and field applications to 
conclude that Landsat digital data is sufficiently accurate to define the 
land cover distributions required as inputs to regional planning and re- 
source allocation models. Even though computer aided translation of raw 
Landsat data is extremely efficient, the adoption of this new technology by 
counties and other regional governements has been limited. A major problem 
continues to center on the difficulties of merging Landsat derived land 
covers into a geographic information system (GIS) that a regional government 
may have been using for a number of years. Typically, the data stored in 
such a GIS is referenced to USGS quadrangle sheets and/or state plane 
coordinates. 

The paper describes an approach for merging multi-scene Landsat data 
bases into existing geographic information systems having 5-second or smaller 
cells. The approach uses the output from the State of Maryland's UNIVAC 
1180-based Landsat classification program ASTEP (Algorithm Simulation Test 
and Evaluation) developed by NASA. The structure of the technique was 
designed to address the problems that emerged as part of the Landsat classi- 
fication of the 64,000 square mile Chesapeake watershed involving twelve scenes 
that was conducted by the senior author as part of an EPA study. The paper 
describes the removal of overlap among adjacent scenes, the crossreferencing 
of ground control points, and the isolation of the appropriate pixels from the 
Landsat data base for subsequent positioning into a file containing ancillary 
data referenced to a specific USGS minute quadrangle sheet. Examples 
illustrate the clustering of classified Landsat pixels to define the dominant 
land use for each of 8, 100 cells within a series o f quadrangle, sheets 
distributed over the State of Maryland. 

The approach uses a hard copy terminal tied to an ASTEP algorithm through 
telephone lines. A coordinate digitizing board for inputing the position of 
ground control points is also valuable, although manual measurements are 
possible. The approach is quite efficient and should be especially attractive 
for use on regional scale studies. 
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1 . 


INTRODUCTION 


Many states, counties, public utilities and other organizations concerned 
with planning and raanagementon a regional scle have integrated computer-based 
geographic information systems (GIS) into their decision making processes. 
Properly designed and operated systems allow the decision maker to define the 
spatial distribution of current conditions within the area of interest and, in 
an increasing number of cases, conditions in surrounding areas. Of equal 
importance, a good GIS also allows the decision maker to better understand how 
the region evolved to its current state and to interpret trends that indicate 
future conditions. When relatively large areas are involved, the use of GIS 
and computer technologies are pivotal in the development of effective planning 
and management strategies. 

Current and past land cover distributions are key elements in the GIS. 
Unfortunately, these land cover files are often poorly defined or not up-to- 
date because of the times and costs required to assemble the data, interpret 
it, encode it and then enter into the GIS. Professionals concerned with GIS 
have long recognized the potential of digital format data from the Landsat 
series of satellites as a base for maintaining up-to-date land cover files in 
their systems. Although there are many successful applications, Landsat has 
remained a "potential" to the typical GIS user because the data format is less 
than ideal. This is especially true for grid cell based GIS that are referenced 
to USGS or state plane coordinate systems. Individual cells in such systems 
typically run from north-south vectors and may be 10, 91.8 or 4.5 acres or they 
may be 5 seconds in size. Although there are a number of programs designed 
to geometrically correct and reformat Landsat data, the time required to learn 
these systems, the level of effort and often the special equipment required 
limits the widespread application of many of the techniques and, thereby, 

Landsat continues to be unrealized potential . 

2. THE CHESAPEAKE BAY EXPERIMENT 

The development of the land cover distributions of the 64,000 square mile 
Chesapeake Bay watershed can be used to illustrate some of the problems that 
regional planning and management organizations encounter when attempting to 
integrate a Landsat derived data base into their operations. Figure 1 shows 
the outline of the Chesapeake watershed and the geometry of the 12 scenes used. 
The objectives of the Chesapeake Bay Project were: 1) produce a Level I land 

cover classification of the Chesapeake Bay watershed; 2) within agriculture 
land cover, determine tillage practices; and 3) tabulate land cover statistics 
by river subbasins. The land cover statistics were required as input to a 
mathematical model to predict the non-point source pollution loads to the 
Chesapeake Bay. The classification was conducted by the Northern Virginia 
Planning District Commission for the Environmental Protection Agency and used 
the IDIMS (Interactive Digital Image Manipulation System) and GES (Geographic 
Entry System) at NASA's Goddard Space Flight Center. The scenes had the 
known geometric distortions corrected (deskewing, removal of synthetic pixels) 
and procedures were developed to remove the overlap among scenes. The result 
was a properly registered land cover distribution that, through the use of a 
digitizer, was summarized for 63 subwatersheds distributed throughout the basin. 
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CHESAPEAKE BAY BASIN 


FIGURE 1 

Landsat Scene Location within Chesapeake Bay Basin 
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The existence of such a digital data base was, obviously, very attractive 
to organizations within the area that had computer-based GIS. State systems, 
such as Maryland's MAGI (Maryland Automated Geographic Information System) and 
county systems, such as MSDAMP (Multi-Scale Data Analysis Mapping Program) used 
by Montgomery County, Maryland appeard to be the logical recipients of the 
data derived in the Chesapeake study. While it was straight forward to extract 
a relatively large polygon such as a watershed from the Chesapeake Landsat data 
base, MAGI and MSDAMP require the definition of land covers within individual 
cells referenced to USGS or state plane coordinates. The general concept of 
MSDAMP is illustrated in Figure 2. MSDAMP is a series of 90 x 90 five second 
cells referenced to 19 USGS 7% minute quadrangle sheets. The user obtains 
information by entering the name of the quadrangle sheet or sheets of interest 
and then extracts information by defining either polygons or individual cells. 
Maryland's MAGI uses either a 91.8 or 4.54 acre cell. Geographical Information 
Systems of the MSDAMP and MAGI types must have one dominant land cover defined 
for each cell in the data base. A schematic of the definition problem is 
illustrated in Figure 3. The domain of a particular USGS quadrangle sheet must 
be isolated from the Landsat data base and then a specific five second cell 
must become computer retrievable to the staff of the user organization. Image 
processing capabilities are available to all Maryland state and local govern- 
mental organizations through the State's UNIVAC 1100 series of computers 
located at the University and State College campuses. As potential State and 
county users of the Maryland portion of the Chesapeake data base moved toward 
integrating this additional information into their GIS, it became obvious that 
the efforts were not going to be widely successful because the needed software 
did not exist in a form that was compatible with UNIVAC 1100 series computers. 
There was no parallel software that could: rotate the Landsat coordinate 

system; reference the individual cells to USGS coordinates; isolate an array 
of cells defining a USGS 7h minute quadrangle sheet and then resample the 
individual pixels to define a single land cover category for a predefined cell 
size. Without such software, the Chesapeake data base provided an excellent 
source of qualitative information, but remained inaccessable to the day-to-day 
user of the established computer-based geographical information systems 
operating within the State. 

3. OBJECTIVES 

If the Chesapeake Landsat-derived data base and similar future Landsat 
efforts are to be integrated into the existing geographical information 
systems, it is necessary to develop additional software to overcome the 
problems discussed above. To be useable, the additional software has to be 
fully integrated into established computer based approaches that are accessable 
and familiar to the users. Because few of the users in the State of Maryland 
have access to color CRT-based interactive image processing systems, the 
software had to be designed to run on a standard UNIVAC 1108 mainframe computer 
and require no more than a modem-connected hard copy terminal for operation. 
Further, because of severe restrictions placed on core storage during the 
daytime hours, the system had to be designed for minimum core storage 
utilization. With these contraints in mind, system development was undertaken 
to meet the following objectives: 
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Figure 2 

Quad-Sheet Storage Arrangement For 
Montgomery County, Maryland 
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Figure 3 

Isolation of a Single Cell from a Landsat Scene 
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1) Develop interactively an equation which relates a Landsat coordinate 
system to a latitude longitude coordinate system. 

2) Create a transformed data base from Landsat imagery compatible with 
preconfigured Geographical Information Systems. 

3) Enter geographic data from a map surface into ASTEP. 

4. SYSTEM CAPABILITIES 

To meet the objectives listed above, two programs (REGISTER and TRANSFORM) 
were developed. The program REGISTER was designed to input geographic data 
from a map surface and develop a regression relating the Landsat coordinate 
system to a latitude, longitude coordinate system. The program TRANSFORM, 
using the equations developed by the program REGISTER, was designed to create 
a geometrically correct data base compatible with preconfigured geographical 
information systems. 

The program REGISTER serves two functions. First, it outputs to a file 
the longitude and latitude of points digitized from a map surface. Second, 
it provides an equation relating the latitude, longitude of a point to its 
Landsat line and sample coordinate. To complete the first function, a link 
(equation) must be developed relating the position on a map surface to its 
latitude and longitude. The position on the map surface may be input as 
coordinates from a digitizing table or measured manually off the map using 
the upper left corner as the origin. (Note: While manually measuring the 

location of points on a topographic map may be tedious, if it is done carefully 
accuracy on a 1:24000 scale map can be ±40 feet.) The user inputs to the 
program are the longitude and latitude in degrees .minutes .seconds of the upper 
left corner of the map, the size of the map in minutes, the distance between 
"tic-marks" on the map in minutes, and the coordinates of the "tic-marks" from 
the digitizing table or as measured manually by the user. A first order 
polynomial regression equation is then developed relating the coordinate from 
the map surface to its longitude and latitude coordinate. A list of the 
actual and predicted coordinates, as well as the residuals, is produced for 
each of the "tic-marks" are output. The user has the option of removing any 
of the "tic-marks" from the registration if they were incorreclty digitized. 

The user also has the option of changing the regression equation to second or 
third order. (Note: For large scale maps, i.e., 1:24000 there should be no 
need to go to a j;ecqnd or third order -equation.-) - Once-the map-has been- 
registered to the digitizing table, the location of the ground control points 
can be digitized from the map. These points are then stored in a file for use 
in developing the transformation equation. 

The second function that the program REGISTER performs is to allow the 
user to develop an equation relating the latitude, longitude coordinate 
system of the Geographical Information System to the line and sample co- 
ordinate of Landsat. The program reads the file containing the ground control 
points created above and a least- square file is applied to the points to 
develop a simple linear transformation of the form: 


X = C(l) + C (2) Y + C(3)X 
Y = C(9) + C(10)Y + C(11)X 

A A 

where X and Y are the estimated sample and line values in the Landsat co- 
ordinate system, X and Y are the observed values of longitude and latitude 
(digitized coordinates) in the GIS coordinate system, and C(l), C(2)...C(N) 
are the coefficients fo the transformation equation expressed in the form: 


C(l) C (9) 

[YX] = [1YX] C (2) C CIO) 

C(3) CCH) 

An output table is printed that contains the estimated sample and line value, 
the observed sample and line value, and the error (observed-estimated) sample 
and line value for each ground control point. 

Upon examination of the ground control points, the user has the option 
of altering the list of ground control points. The user is prompted: 

DO YOU WISH TO EDIT POINTS? Y/N. If the user responds with an upper case 
Y he is prompted with: ADD (A) DELETE (D) OR EXIT(E)?. If the user wishes 

to delete a point, he responds with an upper case D. (Note: the development 
of the equation is an iterative process, the user may wish to restore a 
ground control point that was previously deleted by responding A.) The user 
is then prompted: INPUT NUMBER(S) TO BE ADDED OR DELETED ZERO (0) TO END. 

The user then would input the number(s) of the ground control point to be 
deleted, 0 indicates there are no more points. The user is then prompted: 

ADD(A) DELETE(D) OR EXIT (E) ? and would respond E. The first order regression 
equation is then recalculated and the output table is again listed. The 
user has the option of editing points and recalculating the first order 
regression equation until he is satisfied that all the ground control point 
residuals have the same order of magnitude. When the prompts to edit points 
are answered N, the user will be prompted with: DO YOU WISH THIS TO BE THE 

HIGHEST ORDER? Y/N. If the response is N, a second order equation is developed 
with the form: 

X = C(l) + C(2)Y + C(3)X + C(4)Y 2 + C(5)X 2 + C(6)XY 

Y = C(9) + C(10)Y + C(11)X + C(12)Y 2 + C(13)X 2 + C(14)XY 

The output table is printed and the user is given the option of editing points 
or developing a third order equation. The third order equation has the form: 

X = C(l) + C(2) Y + C(3)X + C(4)Y 2 + C(5)X 2 + C(6)XY + C(7)Y 3 + C(8)X 3 

Y = C(9) + C (10)Y + C(11)X + C(12)Y 2 + C(13)X 2 + C(14)XY + C(15)Y 3 + C(16)X 3 

When the final transformation equation has been calculated, the equation is 
stored in a disc file for use by the program TRANSFORM. 
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The program TRANSFORM reads the transformation equation developed above 
and prompts the user for information concerning the location and size of the 
transformed area. The program prompts the user with: INPUT LONGITUDE, 

LATITUDE OF UPPER LEFT CORNER OF STUDY AREA IN DD.MMSS. When the user 
responds, he is prompted: INPUT THE SIZE OF TRANSFORMED AREA IN MINUTES 

LONGITUDE, LATITUDE. The user is not constrained to having the size of the 
transformed area the same in both longitude and latitude. The user is then 
prompted for the number of cells in the X and Y directions in the transformed 
area. The program TRANSFORM displays the Landsat sample and line value for the 
four comers of the transformed area and the minimum subset of the original 
data needed to transform the area. (Note: This allows the user to redefine 

the original study area to use the minimum amount of computer storage and CPU 
time.) The user then has the option of stopping the run to subset the original 
data, or continuing the run and creating the transformed area. 

Figure 4 is a schematic representation of the procedure to transform a 
study and form a Landsat line and sample coordinate system to a latitude, 
longitude coordinate system. The output from the program TRANSFORM is a 
file that can be read directly into a geographic information system or re- 
formatted by a program (INASTEP) to be entered back into ASTEP to use its 
statistical and map generating capabilities. 

5. PROCEDURE 

The procedure to transform a study area from a Landsat referenced co- 
ordinate system into a georeferenced coordinate system is as follows: 

1) Ouput lineprint maps of study area 

2) Locate and digitize features that can be found on both lineprint 
maps and topographic maps. 

3) Develop regression equation 

4) Transform the data. 

The first step in transforming the data is to output lineprint maps from 
ASTEP, such as that illustrated in Figure 5, of the study are for use in 
locating features (ground control points) . The lineprint map generation is 
the most critical portion of locating ground control points. A lineprint map 
is limited to displaying one channel of data with a practical limit of 20 
grey levels, therefore, whatever a user can do to combine information from 
more than one MSS channel of data on a lineprint map is important. There are 
many ASTEP output products which are useful in the production of lineprint maps. 
A grey level map (density slice) of channel 7 can provide good land/water inter- 
face detection; it can also be useful in locating bridges, river boundaries, and 
power line clear cuts. A grey level map of channel 5 is useful in finding 
man-made features such as road intersections and industrial parks. 

There are three ASTEP routines that allow the user to output information 
from more than one MSS channel of data. A map from the norm of all four 
channels (brightness map) can often be used to agument the output products 
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listed above. An unsupervised classification of the study area using 
relatively few classes (6-8) provides a quick method of locating forest and 
grass field boundaries. Figure 5, for example, is an example of an un- 
supervised classification using only 6 classes. 

Once the lineprint maps of the study area have been generated, the ground 
control point location can begin. It is important to find ground control 
points that are uniformly distributed throughout and surrounding the study 
area. A general rule of thumb is that in order to have confidence in the co- 
efficients, there should be at least four ground control points for each of the 
coefficients of the regression equation. Therefore, a first order equation 
should have a minimum of 12 ground control points, a second order should 
have a minimum of 24 ground control points, and a third order should have a 
minimum of 32 ground control points. Depending on the size of the study area, 
the land cover, and the topography, it may not be feasible to find as many as 
30 ground control points using lineprint maps. The ground control point can 
be any fixed feature locatable on both the lineprint maps and topographic 
maps. They may include bridges, islands, road intersections, power line clear 
cuts, and small ponds. It is generally easier to locate ground control point 
form images in early spring or late fall when there are no leaves on the trees 
to obscure ground features. There are some features however, that are easier 
to locate in summer scenes (i.e., poweT lines, roads). 

When all the ground control points have been located, two files are created 
for each topographic map in the study area. The first file ocntains the 
digitized coordinates of the "tic-marks" on the topographic map. The second 
file contains the digitized coordinates and Landsat sample and line coordinates 
for each of the ground control points. 

After the ground control points have been digitized, the process of 
developing the regression equation can begin. Tables I - V are examples from 
a program runstream which illustrates the process of developing a regression 
equation. For the sake of simplicity, only control points from one topo- 
graphic map will be used. The user reponses are underlined and comments are 
in brackets. 
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TABLE I 


Initial Output Used to Verify Regression Equation Defining 
Coordinates of Points on Quadrangle Sheet 


BXQT RSSL*REGISTER.ABSTS 

INPUT 1 TO DIGITIZED TO DEVLOP REGRESSION, 0 TO QUIT 

1 

ENTER LONGITUDE AND LATITUDE OF THE UPPER LEFT CORNER 
SIZE OF TOPO IN HINUTS(LON,LAT), DISTANCE BETWEEN THE 
REGISTRATION POINTS IN MINUTS(LONrLAT) 

76.3730 39.3000 7,5 7.5 2.5 2.5 

INPUT RER I STRATI ON POINTS 

BADD TRY1. [TRY1. is the file that contains the digitized coordinates of 
the "tic-marks"] 

vj j J w w 


# X/ Y/ 

1 7B. 37300 39.30001 76.37300 

2 76.35001 39.30000 76.35000 

3 76.32300 39.30000 76.32300 

4 76.30001 39.30000 76.30000 

5 76.37299 39.27299 76.37300 

6 76.34599 39.27299 76.35000 

7 76.32300 39.27300 76.32300 

B 76.30000 39.27300 76.30000 

9 76.37300 39.24600 76.37300 

10 76.34599 39.25000 76.35000 

11 76.32299 39.25000 76.32300 

12 76.29599 39.25000 76.30000 

13 76.37302 39.22300 76.37300 

14 76.35000 39.22299 76.35000 

15 76.32301 39.22301 76.32300 

16. 76.29600 39.22299 c 3|^ 0000 

276326 . 42283975470368000000 
.01966168255291389490 

ERROR SQ ; 25:46G6423 Sii gjj}^3}}jjj 3749984 

SUM ERR X s .023436 SUM ERR Y = .015625 

DO YOU WISH TO EDIT POINTS? Y/N 


1 

2 

3. 


Y EX EY 

39.30000 -.0000004 -.0000084 

39.30000 -.0000082 -.0000002 

39.30000 -.0000059 .0000004 

39.30000 -.0000086 -.0000012 

39.27300 .0000121 .0000064 

39.27300 .0000066 .0000070 

39.27300 .0000016 .0000035 

39.27300 -.0000012 .0000000 

39.25000 .0000039 .0000018 

39.25000 .0000066 -.0000018 

39.25000 .0000164 .0000006 

39.25000 .0000113 -.0000010 

39.22300 -.0000219 -.0000051 

39.22300 -.0000016 .0000053 

39.22300 -.0000094 -.0000141 

39.22300 .0000082 

141730.89427985332736000000 
19.76698345728208928000 
- . 00409724894996088550 


DO YOU WISH THIS TO BE THE HIGHEST ORDER? Y/N 

_Y_ 

INPUT GROUND CONTROL POINTS * THEN DIGITIZED COORDINATES AND LABEL 

BADD TRY. [TRY. contains the digitized coordinates of the ground control 
g^QP points and their line and sample coordinates.] 


INPUT 1 TO DIGITIZE, 2 TO DEVELOP REGRESSION, 0 TO QUIT 
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Table I is a list of the predicted and actual longitude and latitude 
of the "tic-marks” on the topographic map, as well as the errors (actual- 

predicted) . 


where : 


XI = predicted longitude of "tic-mark 
Y1 = predicted latitude of "tic-mark" 
X = actual longitude of "tic-mark" 

Y = actual latitude of "tick-mark" 


EX = X - XI 
EY = Y - Y1 

In this example all the errors are less than 0.2 of a second . (approximately 
16’ at this latitude) so there was no need to edit points ^crease te 
rvrvw of the regression equation. After the user replies Y to the prompt 
"DO YOU WISH THIS TO BE THE HIGHEST ORDER?" the user is prompted for the ground 
control points and digitized coordinates. 

The program uses the equations generated in Table I to convert the 
digitUed’coordinates of the ground control points to their correspond! 
i on oi tnde latitude and stores the results for later use. The process is 
repeated for each topographic map in the study area After all map s have een 
digitized, the user responds "2" to the prompt "INPUT 1 TO DIGITIZE, 2 TO 
DEVELOP REGRESSION, 0 TO QUIT". 

TABLE II 

Output Used to Verify Regression Equation 


f X/ Y / X 

1 1598.32462 125.30842 1600.00000 

2 1700.71982 121.36528 1707.00000 

3 1752.95477 102.56207 1767.00000 

4 1673.43864 162.95376 1670.00000 

5 1766.56570 142.39142 1772.00000 

6 1800.19225 214.15857 1794.00000 

7 1749.11740 170.16807 1749.00000 

8 1811.80373 223.93128 1804.00000 

9 1793.64261 169.85357 1795.00000 

10 1693.76706 247.05843 1674.00000 

11 1730.37199 283.31932 1775.00000 

12 1799.00040 254.53150 1783.00000 

13 1764.99077 209.58345 1758.00000 

14 1718.71126 179.44728 1714.00000 

15 1678.39890 195.36758 1670.00000 

# X C0EFF 

1 13354 5:i??OT?r“ 

3 -.36905888861332414208 

ERROR SQ = 3226.59675617842002880000 

SUM ERR X= -000092 SUM ERR Y = .000011 


126.00000 

121.00000 

103.00000 

164.00000 

143.00000 

212.00000 

170.00000 

223.00000 

169.00000 

245.00000 

290.00000 

252.00000 

209.00000 

180.00000 
195.00000 

38163. 


EX 

1.6753845 
6.2801819 
14.0452271 
-3.4386444 
5.4342957 
-6.1922455 
-.1174011 
-7.8037262 
1.3573914 
-19.7670593 
44.6280060 
-16.0003967 
-6.9907684 
-4.7112579 
-8.3988953 
Y COEFF 


EY 

.6915751 

-.3652592 

.4379263 

1.0462418 

.6085815 

-2.1585655 

-.1680679 

-.9312840 

-.8535690 

-2.0584335 

6.6806755 

-2.5314999 

-.5834541 

.5527229 

-.3675785 


.06160473314834291592 
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The program reads the file containing the latitude, longitude and line 
sample values for each ground control point and develops a regression equation 
relating latitude, longitude to line and sampple. Table II is a list of the 
actual line and sample^ predicted line and sample and the errors for each of 
the ground control points. Where: 

XI = predict sample value 

Y1 = predict line value 

X = actual sample value 

Y = actual line value 

EX = X = XI 

EY = Y - XI 

The user is prompted "DO YOU WISH TO EDIT POINTS? Y/N". Deciding which 
point (s) to remove from the regression equation is somewhat of an art. A good 
rule-of-thumb would be to remove any point whose errors are significantly 
different from the rest (i.e., point 11). Being an iterative process, the user 
can delete points to see the effects and later add them if he wishes. In 
this example, the point "11" is deleted. 

TABLE III 


Output Used to Verify Regression Equation Without Point #11 

DO YOU WISH TO EDIT POINTS'? Y/N 


ADD (A) DELETE(D) OR EXIT(E) ? 


INPUT NUNBER(S) TO BE ADDED OR DELETED ZERO ( 0 ) TO END 


11 


ADD (A) DELETE (D) 

OR EXIT (E) 

7 




* 

X/ 

Y/ 

X 

Y 

EX 

EY 

1 

1597.25696 

125.22393 

1600.00000 

126.00000 

2.7430420 

.7760677 

2 

1706.68803 

122.18216 

1707.00000 

121.00000 

.3119659 

-1.1921606 

3 

1765.75558 

104.30714 

1767.00000 

103.00000 

1.2444153 

-1.3071394 

4 

1669.70222 

162.44365 

1670.00000 

164.00000 

.2977753 

1.5563526 

5 

1772.50732 

143.16842 

1772.00000 

143.00000 

-.5073242 

-.1684227 

6 

1794.33221 

213.26280 

1794.00000 

212.00000 

-.3322144 

-1.2628002 

7 

1748.62469 

170.05817 

1749.00000 

170.00000 

.3753052 

-.0581665 

8 

1804.76753 

222.86302 

1804.00000 

223.00000 

-.7675323 

.1369820 

9 

1795.93893 

170.10218 

1795.00000 

169.00000 

-.9389343 

-1.1021824 

10 

1675.03110 

244.43800 

1674.00000 

245.00000 

-1.0310974 

.5619984 

12 

1785.26913 

252.53663 

1783.00000 

252.00000 

-2.2691345 

-.5366306 

13 

1757.85748 

208.53445 

1758.00000 

209.00000 

.1425171 

.4655514 

14 

1714.56311 

178.84788 

1714.00000 

180.00000 

-.5631104 

1.1521187 

15 

1668.70555 

194.02156 

1670.00000 

195.00000 

1.2944489 

.9784431 

# 


X 

: C0EFF 


Y COEFF 



1 134280. 

80787729471872000000 

38173.89989193249484800000 


2 -.15026574780891621632 

3 -.40357207301249786560 

ERROR SQ = 31.58291904046200204800 

SUM ERR X= .000122 SUM ERR Y = .000011 


-.37825928840648126272 

.05704062922797348256 
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Table III lists the output for the regression equation developed without 
ground control point "11". There are no ground control points having errors 
significantly different from the rest, so the user responds N to the prompt 
"DO YOU WISH TO EDIT POINTS? Y/N". The user is then prompted with "DO YOU 
WISH THIS TO BE THE HIGHEST ORDER? Y/N". If the user is not satisfied with 
the size of the errors, he will respond "N" and a second order equation will 
be developed, as illustrated in Table IV. 

TABLE IV 

Output to Verify Regression Equation Using Second Order Equation 


DO 

1 

YOU NISH TO EDIT POINTS? 

1 Y/N 




1 

DO 

| 

YOU WISH THIS 

TO BE THE 

HIGHEST ORDER? Y/N 



1 

# 

X/ 

Y/ 

X 

Y 

EX 

EY 

1 

1599.11740 

126.62065 

1600.00000 

126.00000 

.8825989 

-.6206522 

2 

1707.99907 

122.20280 

1707.00000 

121.00000 

-.9990692 

-1.2027960 

3 

1767.02563 

103.49714 

1767.00000 

103.00000 

-.0256348 

-.4971447 

i 

14?S:IS§il 

113 : 11 ??? 

1998:88888 

mmi 


:!}?§??? 

6 

1793.28535 

212.64375 

1794.00000 

212.00000 

.7146454 

-.6437511 

7 

1748.69167 

169.81061 

1749.00000 

170.00000 

.3083344 

.1893864 

8 

1803.47786 

222.15224 

1804.00000 

223.00000 

.5221405 

.8477650 

9 

1795.66638 

169.25774 

1795.00000 

169.00000 

-.6663816 

-.2577400 

10 

1674.15500 

245.56256 

1674.00000 

245.00000 

-.1549988 

-.5625591 

12 

1783.61646 

252.21816 

1763.00000 

252.00000 

-.6164551 

-.2181587 

13 

1757.14552 

208.36650 

1758.00000 

209.00000 

.8544769 

.6335011 

14 

1714.68178 

179.08262 

1714.00000 

180.00000 

-.6917760 

.9173832 

15 

1668.79240 

194.93917 

1670.00000 

195.00000 

1.2075958 

.0608253 

I 

1 


X COEFF 


Y COEFF 



1 55196.07290533184998400000' 8554.40744912996984480000 

2 -.73650654382436186816 .15635408307161924240 

3 .46676819108142608384 -.00889545607060426843 

4 .00000245849135216414 -.00000080964223510263 

5 -.00000147942547933447 .00000041454940314387 

6 -.00000037657992861539 -.00000110603672397266 

ERROR SQ * 13.09050019155256448000 

8UM ERR X s .008270 SUM ERR Y = .005571 

DO YOU WISH TO EDIT POINTS? Y/N 


DO YOU NISH THIS TO 8E THE HIGHEST ORDER? Y/N 


END PROGRAM REGISTER. 
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The output shown in Table IV is again listed for the second order 
equation and the user has the option to edit points or go on to a third order 
equation. Once the user responds "N" to each question, the program creates 
a file with the coefficients of the regression equation. 

After the regression equation has been developed, the user can transform 
the study area. The program TRANSFORM using the regression equation developed 
above, prompts the user for location and format of the transformed area. The 
program reads the raw Landsat data from a file created by ASTEP, transforms 
the data and outputs the transformed data in a format compatible with various 
geographic information systems. The user will be prompted for the longitude, 
latitude of the upper left corner of the study area, and the size of study 
area in minutes. The user is not required to have the study area correspond 
to one topographic map, and the study area can have different dimensions in 
the latitude and longitude direction. Table V is an example run of the program 
TRANSFORM; the study area is the Towson, MD quadrangle. (See Figure 6) The out- 
put file is to have the data stored in 5 second cells. 

TABLE V 

Example Row for Towson, MD 


INPUT LON.LAT OF UPPER LEFT CORNER OF AREA IN D.MS 
76.3730 39.3 000 

INPUT SIZE OF STUDY AREA IN MINUTS LON.LAT 
7.5 7.5 

INPUT NUMBER OF CELLS LON.LAT IN TRANSFORMED AREA 
90 90 

LANDSAT COORDINATES OF T0P0 SHEET 

1590 .. 121 .**#**##*## 1770 . . 93 . 

I I 

I I 

1 I 

I - - - - I - 

1655 ., 292 .*****»»*** 1934 . » 264 . 

IF YOU NISH TO SUBSET STARTING LINE STARTING SAMPLE# LINES # SAMPLES 

66 1594 209 256 

DO YOU NISH TO SUBSET ? Y/N 
Y 

DO YOU NISH TO QUIT ? Y/N 

JL 

END PROGRAM TRANSFORM 
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FIGURE 6 

USGS 7,5' Topographic Map Towson, MD 
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Figure 5 is the raw Landsat data used in the example in Table V, and Figure 7 
is a map of the output file. 

The user, by changing the number of cells in the output transformed area 
can produce output products at a given scale (i.e., 1:24000). If in the example 
above, had changed the number of cells in the transformed area from 90 x 90 
to 212 x 272, a map of the output file would have a scale of 1:24000. Figure 
8 is an example of such a product. 

6. CONCLUSION 

Many current or potential users of digital format remotely sensed imagery 
are restricted to the use of a remote lineprinter type terminal that accesses 
processing software on a general purpose, mainframe computer. The software 
described in the present paper was designed to provide this group of users with 
some of the interactive geometric corrections and data manipulation capabilities 
found on dedicated, color CRT-based image processing systems such as IDIMS. The 
system developed is compatible with ASTEP input/output reoutines and the UNIVAC 
1100 series core limitations. It requires only a typewriter type terminal and 
is, therefore, available to Maryland State and local government users. 

The interactive editing capabilities allow the user to produce a ±1 pixel 
registration accuracy between an image and map referenced position. Flexible 
output formate routines allow interfacing with preconfigured geographical 
information systems. With minor modifications, the system can easily be adapted 
to other geographical formats (i.e., state plane, UTM) and other sensors 
(i.e., RBV) . The resulting transformed data bases can be re-entered into the 
ASTEP program to allow the user access to ASTEP capabilities such as scaled 
map production and statistical tabulations. 
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CLASSIFIED IMAGE 


TOWSON MD. QUADRANGLE 
GEO-CORRECTED 90x90 5 SECOND CELLS 


Figure 7 
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FIGURE 8 

Quad-Centered Transformation Scale 1:24000 


181 



